Bilayer superfluidity of fermionic polar molecules: many body effects 



O 

(N 

o 

Q 

(N 



C/2 

^— > 



I 

c 

o 
o 



> 

On 
00 

in 
in 



X 



M.A. Baranov, 1 - 2 ' 3 A. Micheli, 1 ' 2 S. Ronen, 1 ' 2 and P. Zoller 1 - 2 

institute for Theoretical Physics, University of Innsbruck, A-6020 Innsbruck, Austria 
2 Institute for Quantum Optics and Quantum Information of the Austrian Academy of Sciences, A-6020 Innsbruck, Austria 
1 1 RRC "Kurchatov Institute", Kurchatov Square 1, 123182 Moscow, Russia 
(Dated: December 30, 2010) 

We study the BCS superfluid transition in a single-component fermionic gas of dipolar particles 
loaded in a tight bilayer trap, with the electric dipole moments polarized perpendicular to the layers. 
Based on the detailed analysis of the interlayer scattering, we calculate the critical temperature of 
the interlayer superfluid pairing transition when the layer separation is both smaller (dilute regime) 
and of the order or larger (dense regime) than the mean interparticle separation in each layer. 
Our calculations go beyond the standard BCS approach and include the many-body contributions 
resulting in the mass renormalization, as well as additional contributions to the pairing interaction. 
We find that the many-body effects have a pronounced effect on the critical temperature, and can 
either decrease (in the very dilute limit) or increase (in the dense and moderately dilute limits) the 
transition temperature as compared to the BCS approach. 
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I. INTRODUCTION 

Recent experiments have prepared quantum degener- 
ate gases of homonuclear and heteronuclear molecules 
in electronic and vibrational ground states [H-Q. Het- 
eronuclear molecules, in particular, have large electric 
dipole moments associated with the rotational excita- 
tions, The new feature of polar molecular gases is thus 
the strong, anisotropic dipole-dipole interactions between 
the molecules, which can be controlled with external elec- 
tric fields [7Hl0lj . When DC electric fields are applied to 
polarize molecules, a major obstacle is given by the in- 
crease of inelastic collision rates corresponding to chem- 
ical reactions between the molecules, as in the case of 
KRb in the recent JILA experiments. However, these 
can be strongly suppressed, and thus the gas stabilized 
by tightly confining the molecules in a single plane of 
a quasi-2D geometry [ll|. This relies on the fact that 
for induced electric dipole moments perpendicular to the 
plane of confinement, the dipolar forces will be repulsive 
in-plane, thus suppressing short distance inelastic colli- 
sions. Such a 2D confinement can be achieved by loading 
the gas of polar molecules into a ID optical lattice. This 
leads naturally to a setup of a multilayer polar gas where, 
however, forces between dipoles in different layers can be 
attractive, and the collapse is prevented by a sufficiently 
high optical potential barrier. For a bilayer system this 
attraction can lead to the formation of bound pairs of 
polar molecules, reminiscent of bilayer excitons, and in 
a multilayer configuration to the formation of strings of 
molecules. In particular, in a gas of single component 
fermions loaded in a tight bilayer trap, as realized with 
KRb, this will give rise to an s-wave BCS superfluid tran- 
sition [l2| (for p-wave pairing in a monolayer of polar 
molecules see Refs. [l]| and • It is the purpose of 
this work to study this bilayer BCS superfluid transition 
in some detail, in particular we go beyond Ref. [l2[ with 
emphasis on the inclusion of many body effects. 



In the bilayer BCS superfluid the single species polar 
molecules in the two layers provides the system with a 
two-component character, where two species are particles 
on different layers coupled by the long-range dipole in- 
teraction, allowing fermions from different layers interact 
in the s-wave channel that is dominant at low energies. 
This interlayer interaction results in very peculiar prop- 
erties in both a two-body system: the existence of in- 
terlayer bound states lla| -|18| and various regimes of the 
interlayer scattering 1181. and in a many-body system: 
interlayer BCS pairing[12| and BCS-BEC crossover fl2j |. 
[l9j . Based on a detailed analysis of various regimes of 
interlayer and intralayer scattering, we extend the anal- 
ysis of the interlayer superfluid pairing beyond the stan- 
dard BCS approach [12:] by including many-body contri- 
butions resulting in the mass renormalization, as well as 
in additional contributions to the pairing interaction. We 
perform the calculation of the critical temperature in the 
regime of a weak interlayer coupling for the cases when 
the layer separation is both smaller (dilute regime) and 
of the order or larger (dense regime) than the mean in- 
terparticle separation in each layer. As found, the many- 
body effects have a pronounced effect on the critical tem- 
perature, and could either decrease (in the very dilute 
limit) or increases (in the dense and moderately dilute 
limits) the transition temperature as compared to the 
BCS approach. 

The paper is organized as follows. Sec. II gives an 
overview of the problem and identifies the relevant pa- 
rameters and parameter regimes. In Sec. Ill we intro- 
duce the model describing bilayer pairing. Sec. IV dis- 
cusses two-particle bound states and scattering proper- 
ties for the interlayer problem. Sec. V provides a theo- 
retical treatment of many body effects in interlayer BCS 
pairing. Results for the critical temperature in the dilute 
and dense limit are summarized in Sees. VI and VII, 
respectively. 
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II. OVERVIEW OF THE PROBLEM AND 
SUMMARY OF THE RESULTS 

The considered single-component fermionic bilayer 
dipolar system provides an example of a relatively simple 
many-body system in which an entire range of nontriv- 
ial many-body phenomena are solely tied to the dipole- 
dipole interparticle interaction with its unique proper- 
ties: long-range and anisotropy. The long-range charac- 
ter provides an interparticle interaction in single com- 
ponent Fermi gases inside each layer which otherwise 
would remain essentially noninteracting. For the consid- 
ered setup, this intralayer interaction is always repulsive 
and gives rise to the crystalline phase for a large den- 
sity of particles. More important, the long-range dipole- 
dipole interaction couples particles from different layers 
in a very specific form resulting from the anisotropy of 
the interaction: two particle from different layers attract 
each other at short, and repel each other at large dis- 
tances, respectively. The potential well at short distances 
is strong enough to support at least one bound state for 
any strength of inter layer coupling. For a weak coupling 
between layers, the bound state is extremely shallow and 
has an exponentially large size. However, in the inter- 
mediate and strong coupling cases the size of the deepest 
bound state becomes comparable with the interlayer sep- 
aration. In the fermionic many-body system this behav- 
ior of the interlayer interaction leads to a BCS state with 
interlayer Cooper pairs in the weak (interlayer) coupling 
regime when the size of the bound state is larger than the 
interparticle separation (in other words, the Fermi energy 
is larger than the binding energy) . With increasing inter- 
layer coupling, this BCS state smoothly transforms into 
a BEC state of tightly bound interlayer molecules when 
the interparticle separation is larger than the size of the 
bound state. Of course, the BEC regime and BEC-BCS 
crossover are only possible when the mean interparticle 
separation in each layer is larger than the distance be- 
tween the layers. 

In this paper we focus on the regime of a weak intra- 
and interlayer interactions, that allows the usage of con- 
trollable calculations on the basis of the perturbation the- 
ory, and consider in details the formation of the interlayer 
BCS state. Before entering the technical derivation, it 
seems worthwhile to briefly identify the relevant param- 
eters and parameter regimes for the bilayer many body 
system of single species fermionic dipoles. In addition, 
we will point out, where the relevant results and discus- 
sions for these parameter regimes can be found in later 
parts of the paper. 

It follows from the previous discussion that the con- 
sidered system is characterized by three characteristic 
lengths: the dipolar length ad = rnd 2 /ti 2 , where m is 
the mass of dipolar particles with the (induced) dipole 
moment d, the interlayer separation I, and the mean 
interparticle separation inside each layer ~ fc^ 1 with 
kp = V \-nn being the Fermi wave vector for a 2D single- 
component fermionic gas with the density n. Therefore, 



the physics of the system is completely determined by 
two dimensionless parameters which are independent ra- 
tios of the above lengths. 

The first parameter g is the ratio of the dipolar length 
and the interlayer separation, g = ad/ 1, and is a measure 
of the interlayer interaction strength relevant for pairing. 
In experiments with polar molecules, the values of the 
dipolar length ad is of the order of 10 2 -j- 10 4 nm: for 
a 40 K 87 Rb with currently available d 0.3 D one has 
ad ~ 170 nm (with ad ~ 600 nm for the maximum value 
d « 0.566 D), and for 6 Li 133 Cs with the tunable dipole 
moment from d — 0.35 D to d = 1.3 D in an external 
electric field ~ 1 kV/cm the value of ad varies from ad ~ 
260 nm to ad ~ 3500 nm. For the interlayer separation 
I = 500 nm these values of ad corresponds to g < 10. 

The second parameter kpl measures the interlayer sep- 
aration in units of the mean interparticle distance in each 
layer. This parameter can also be both smaller (dilute 
regime) and of the order or larger (dense regime) than 
unity for densities n = 10 6 -j- 10 9 cm -2 (for, example, for 
I = 500 nm one has kpl = 1 for n « 3 • 10 7 cm~ 2 ). 

The two parameters g and kpl determine the regime 
of interlayer scattering at typical energies of particles 
(~ Fermi energy £p = K k F /2m), and their product, 
gkpl = adkp, being the ratio of the dipolar length 
and the mean interparticle separation inside each layer, 
controls the perturbative expansion in the system and, 
therefore, many-body effects. The existence of different 
regimes of the interlayer scattering originates from two 
characteristic features of the interlayer dipole-dipole in- 
teraction, as discussed in the context of Eq. ([JJ below: 
the presence of the typical range ~ / beyond which the in- 
teraction is attactive, and of the long-range dipole-dipole 
repulsive tail. As a result, the Fourier component of the 
interaction (see Eq. flSJ) below), decays exponentially for 
large momentum k 3> l~ l , while it is proportional to 
k for k -C I -1 , and, hence, vanishes for k = 0. This 
leads to three different regimes of scattering and, there- 
fore, of the BCS pairing, depending on the relation be- 
tween g and kpl: regime a when g < kpl < 1, regime 
b when exp(— 1/g 2 ) <C kpl < g < 1, and regime c when 
cxp(— l/g 2 ) < kpl -C g < 1. The exponential factor 
in the last two formulae is related to the size of an ex- 
tremely shallow (in the limit of small g) interlayer bound 
state (compare see Eq. (TTTj) below). 



A. Regime a: g < kpl < 1 

In this regime g < kpl < 1, and the scattering is dom- 
inated by the first Born approximation (see Eq. (|18[) and 
([2"5)l for the s-wave scattering amplitude for kpl -C 1 and 
^fI J; 1) respectively). The critical temperature in this 
case will be given below in Eqs. ([57)) and ([75]) for the di- 
lute, kpl <C 1, and the dense, kpl < 1 cases, respectively. 
The ratio of the critical temperature to the Fermi energy 
(chemical potential), T c /sp, can reach in this regime val- 
ues of the order of 0.1, see Fig.[7J making an experimental 
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realization of the interlayer BCS pairing very promising. 
Note that the maximum of the ration T c /sp corresponds 
to « 0.5. 



B. Regime b: exp(-l/g 2 ) < k F l < g < 1 

In the regime b, exp(— 1/g 2 ) <C kpl < g < 1, the inter- 
layer scattering is dominated by the second order Born 
contribution, see Eq. (|19p . and the critical temperature 
will be given below in Eq. (|62|) . 

C. Regime c: exp(— 1/g 2 ) < HfI < j < 1 

Finally, in the regime c, exp(— 1/g 2 ) < kjrl <C g < 1, 
we recover the universal behavior for a two-dimensional 
low-energy scattering, see Eq. (f20|) . with the typical 
inverse-logarithmic dependence of the s-wave scattering 
amplitude on the energy. The critical temperature is pro- 
vided by Eq. (|63[) and coincide with the critical temper- 
ature in a two component gas with a short-range inter- 
action with the same Fermi energy and the bound state 
energy. Note that the values of the critical temperature 
in the regimes b and c are much smaller (T c //i < 10~ 3 ) 
than in the regime a. This makes an experimental re- 
alization of the interlayer pairing in these regimes very 
challenging. 




FIG. 1: The setup of the dipolar bilayer system: two lay- 
ers with the thickness lo of a one-demensional optical lattice 
formed by two counterpropagating laser waves with wave vec- 
tors \ll and — kz,, are filled with dipoles oriented perpendic- 
ular to the layers. The interlayer distance I is Tr/fci. An 
interlayer Cooper pair/molecule is schematically indicated by 
the dashed oval. 



III. THE MODEL 



We consider a system of single-component polarized 
fcrmionic dipolar particles (harmonically) confined in two 
infinite quasi-two-dimensional layers separated by a dis- 
tance I, which is much larger than the confinement length 
lo in each layer, I ^> lo- We assume that each layer has 
the same density n of dipolar particles with mass m and 
dipolar moment d polarized along the z-axis, which is 
perpendicular to the layers (see. Fig. [I}. The Hamilto- 
nian of the system reads 

H=J2 J ^t (r) |_^. A +ima;X- / i'}^(r) 

+ \ W rfrdr'^t ( r )^t (r ^ (r _ r')i, (r')i> a (r), (1) 

where a = ± is the layer index, z± = z ± 1/2, ipa{*) 
with r = (p, z) is the field operator for fermionic dipolar 
particles (p = xe x + ye y ) on the corresponding layer a, 
A = A p + d 2 /d z is the Laplace operator, uj z is the con- 
fining frequency in each layer such that Iq = ^h/mui z , 
and fi' is the chemical potential. The last term with 



, s d 2 ( z 

V(r) =1-3- 



describes the intra- (a — /3) and interlayer (a =/= j3) 
dipole-dipole interparticle interactions. Assuming a 
strong confinement, ftw z 3> ^',T, where T is the tem- 
perature, we can write 



(V^o) 1/2 



and, therefore, reduce the Hamiltonian (fTJ) to 
H 2D =Y, j d P^(p) {~^ A p ~ Mp) (2) 

+ / dpdp'^i{p)^{p')v a0 {p- P ')Mp')Mp), 



for a two-component fermionic field ?p a (p), a = ±, with 
shifted chemical potential fi = // — Tuu z /2. The intra- 
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component (intralayer) interaction is 
V aa (p) = J dzdz'V{r - Y')^l{z a )4>l{z' a ) 



o J o 



■ exp 



' 1 1 



exp 



4Z 2 



2 + f)K (|)-| Kl (| 



(3) 



where K„(z) is the modified Bessel Functions, and the 
intercomponent (interlayer) one 



V + -(p) = V-+(p) = V 2D (p)^d 2 



2F 



( j0 2 + ;2)5/2- 

The intralayer interaction V++ is purely repulsive 



(4) 



V++(p) 



d 2 /p 3 for p > Z 

{y/2fr){JP/l%)Mlo/p) for p<f(J 



As a result, if the density n is not too large (such that par- 
ticles in each layer are in a gas phase [2U |. (HI), ^ leads 
only to Fermi liquid renormalizations of the parameters 
of the Hamiltonian @ (effective mass, for example). The 
corresponding Fourier transform reads 



6 to 



(5) 



where 

n+(k) 



-27rd 2 fcexp(fc 2 lo/2)[l - erf (klo/Vty] 



with erf (z) = (2/y/n) J? dse s being the error function. 
In the considered limit klg <C 1, one simply has 



-2nd 2 k 



2ttH 2 



-gkl. 



(6) 



The effect of the interlayer interaction V2d(p) (see Fig. 
[5]) is more interesting. A peculiar property of V 2 d(p) is 



-ql 



(7) 



(8) 



dpV 2D (p) -0. 

This means that its Fourier transform 

V 2D (cO = [ dpV 2D (p)e-^P = -—gqle 
J m 

vanishes for small q, 

V 2D (<1^ 0) pa gql -> 0. 



Hence, for sufficiently small energies, the interparticle 
scattering will be dominated by higher order contribu- 
tions in the Born expansion. Note also that the linear 




FIG. 2: The interlayer potential V+-(p) = V2d{p) and its 
Fourier transform ^20(5). 



dependence of V 2 o (q) on q for ql <C 1 is the consequence 
of the long-range power decay of V 2 d (p) for large p (the 
so-called anomalous contribution to scattering, see |23|). 
Another consequence of Eq. ([7]) is that the interlayer 
bound state, which always exists in the potential V 2 d(p) 
[l5j]-[l8j], has extremely low binding energy for g <C 1. 

From the point of view of many-body physics, however, 
the crucial observation is that V 2 d(q) is negative for all 
q, signalling the possibility of the interlayer pairing in the 
form of BCS pairs when the size of the interlayer bound 
state Rb is much large that the interparticle separation, 
Rb ^> nT 1 / 2 , or in the form of interlayer dimers for Rb < 
nT x l 2 with some crossover in between (analogous to the 
BEC-BCS crossover in 2D and 3D for two-component 
Fermi gases with short-range interactions). 

The Hamiltonian ([2} is characterized by two parame- 
ters: g = md 2 /h 2 l, which is the ratio of the dipolar length 
ad = md 2 /h 2 and the interlayer spacing I, and fep^where 
kp = (27m) 1 / 2 is the Fermi wave vector {pp — hkp is 
the Fermi momentum), which is the ratio of the inter- 
layer spacing I and the average interparticle separation 
in each layer. Note that the quantity gkpl = aakp mea- 
sures the strength of the intralayer dipole-dipole interac- 
d 2 kp to the mean kinetic energy 



tion energy d 2 r?l 2 

of particles ~ ep — h 2 kp/2m. In this paper we consider 
weakly interacting gas of dipolar particles with gkpl < 1 
in the dilute kpl < 1 and dense kpl > 1 regimes. 



IV. THE INTERLAYER TWO-BODY 
PROBLEM: BOUND STATES AND 
SCATTERING PROPERTIES 

Let us first discuss the bound state and the scattering 
of two particles interacting with the potential V 2 d{p) - 
interlayer two-body problem (this problem was also ad- 
dressed in Ref. 



18|, see, however, discussion below). For 



this purpose we have to solve the 2D Schodinger equation 
for the relative motion wave function tp(p) 



2m r 



A P + V 2D (p)}ip(p)=ETP(p), 



(9) 
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where m r — m/2 is the reduced mass and the function 
ip(p) is regular for p — > 0. For the bound state solution 
with E = —E b < 0, the wave function ip(p) should decay 
exponentially for large p, while for the scattering wave 
function ip^ip) with E = h 2 k 2 /m > 0, the boundary 
condition for large p reads 



where 7 ~ 0.5772 is the Euler constant, and with the size 



4 +] 



(p) w exp(ikp) - 



^/Snikp 



exp(ifcp), 



where fki'p) is the scattering amplitude and ip is the az- 
imuthal angle, hp = kpcos(ip). (Our definition of the 
scattering amplitude correspond to that of Ref. [23j j . 
which is differs by a factor of — -\/87rfc from the defini- 
tion of Ref. H3.) 



Bound state 



R b = ^h 2 /mE b ~ I exp(4/ ff 2 ) > I. 



Note that the expression ([TT|) for the binding energy co- 
incides with the corresponding the one given in Ref. (l8j 
only to the leading order (~ 1/g 2 )- This is because the 
next order terms in the exponent (~ g^ 1 and ~ g°) are 
determined by the terms of the third and fourth order in 
g in the scattering amplitude (or in the wave function), 
respectively (see Appendix A). In Ref. [IH however only 
terms up to second order were taken into account and, 
therefore, only the leading term is correct, see Figs. [S] 
andEB 

B. Scattering 



Writing the wave function for the relative motion of 
two particles in the form 

■0(p) = Xm,(p)exp(irn z (^), 

where m z is the magnetic quantum number, we ob- 
tain the following equation for the radial wave function 
Xm z (p) of the bound state with the binding energy E b : 



d 2 ld_ 

dp 2 p dp 



mj _ E b + V 2D (p) 

r,2 



h 2 / 



III 



Xmz (p)=0, (10) 



and the function Xm z (p) should be regular for p — > and 
decays exponentially for large r. 

It is sufficient for our purposes to consider the az- 
imuthally symmetric case m z = 0, for which we con- 
sider two limiting cases g ^> 1 and g <C 1. In the 
first case, the potential gv(r) supports several (~ g 1 ^ 2 ) 
bound states, and the lowest bound state has the bind- 
ing energy e b w 2(1 — \J"i/g) and size r b ~ (12g) _1//4 
(E b = {h 2 /ml 2 )2g{\ - yffig) and R b ~ ^Ug)- 1 ^ in 
normal units). In the opposite limit j< 1, there is only 
one shallow bound state (the existence of this bound state 
was proven in Ref. fl5| ) with the binding energy (see de- 
tails of the derivation in Appendix A) 



E h 



Ah 2 
ml 2 



exp 



128 
15g 



2521 
150" 



27 + 0(5) 



(11) 



For the analysis of scattering it is convenient to intro- 
duce the vertex function T(E, k, k'), where the arguments 
E, k,and k' are independent of each other. This function 
satisfies the following integral equation 



r(£,k,kO = F 2D (k-k') + J j^v 2D (k-q) 



1 



E - H 2 q 2 /m + iO 



r(£,q,k'). (12) 



The off-shell scattering amplitude 



/k(k') = 



W- 



dpexp(-zk'p)^ 2D (p)^ +) (p) 



with k k' can be obtained from (m/h 2 )T(E, k, k') 
by putting E — h 2 k 2 /m, and the scattering amplitude 
fk (</?), where if is the angle between k and k', corresponds 
to (m/h 2 )T(E,k,k') with E = H 2 k 2 /m = h 2 k' 2 /m. 



The iteration of Eq. ([12]) up to the fourth order term 
in V2£>(q) reads 



dq V 2D {k- q) V 2D (q - k') f dqidq 2 ^(k - qi)Vfo(qi - q2)F 2 r>(q 2 - k') 



(2tt) 4 (E - h 2 ql/m + i0){E - h 2 q 2 /m + iO) 

dqidqgdqg V 2 p{k - qi)y 2 g(qi - q2)V2p(q2 - <te)V2p(q3 - k') 
(2tt) 6 (E - h 2 q 2 /m + i0)(E - h 2 q 2 /m + i0)(E - h 2 q 2 /m + iO) 

r (1) (E, k, k') + r (2) (E, k, k') + r (3) (e, k, k') + r (4) (e, k, k') + . . . . (13) 
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We now estimate the leading contributions of these 
terms in the small energy limit k ~ k' ~ ^JmE/W < 1/1: 



rW( J B,k,k') = ^(k-k') 



27rn 2 


9 2 


m 


4 ' 




4<? 3 


TO 


15 


2tt?i 2 




m 


32 1 



|k-k'|Z, (14) 
(15) 
(16) 

+ wr]. (17) 



The estimate for rW is trivial, the leading contribu- 
tions to I^ 2 ) and I^ 3 ) come from large q (g> k) and large 
<Zl) 92 3> k) regions, respectively, and the leading con- 
tributions to I^ 4 ) originates from large qi (qi k) and 
<?3 (<Z3 ^ k) but small 52 (92 ~ \JmE/h 2 ). Note that 
the next order terms (except those for I^ 1 )) have relative 
magnitude of the order of (kl) 2 ln(fcZ). 

As already noted in Sec II. the above estimates show 
that there are three different regimes of scattering for 
g < 1 and kl < 1 (dilute weakly interacting regime for a 
many body fermionic system with k ~ k E and E ~ ep ): 
a. g < kl < 1, b. exp(— 1/g 2 ) <C fci < .9 < 1, and c. 
exp(-l/.g 2 ) < fc/ < fif < 1. 

Regime a: The leading contribution to scattering is 
given by the first Born term 



27T^ 2 



valid for g < kl < 1. 

Regime b: The scattering in this case is dominated by 
the second order Born contribution 

r h (E, k, k') « r< 2 > (£, k, k') « - ^ £ (19) 



valid for exp(— 1/3 2 ) <C kl < g < 1. In this case 
the scattering amplitude is momentum and energy in- 
dependent and, hence, is equivalent to a pseudopotential 
V (p) = ~(2nh 2 /m)(g 2 /4)6(p). 

Regime c: In this case the second and the fourth order 
contributions become of the same order and one has to 
sum leading contributions from the entire Born series. 
The result of this summation is 



r c (£,k,k') 



2i:h 2 



to ln(-Ef) /E) + Mr 



(20) 



valid for exp(—l/g 2 ) < kl <C g < 1, where Eb is the en- 
ergy of the bound state from Eq. (jlip . This expression 
recovers the standard energy dependence of the 2D low- 
energy scattering. The scattering amplitude T C (E, k, k') 
has a pole at E = —Eb, as it should be, and the real 
part of T C (E, k, k'), being zero at E = Eb, changes from 
negative to positive values for E > Eb and E < Eb, re- 
spectively. Note that within the lowest order terms, one 



can write a unique expression for the scattering ampli- 
tude for the three regimes in the form (for more details 
see Appendix A) 



T{E, k,k') 



2-nh 2 



3|k-k'| 



HEb/E) 



For later discussion we note that the scattering am- 
plitude has both real and imaginary parts. The relation 
between them can be established on the basis of Eq. (Tl"2"|) 
by considering the imaginary part of both sides of this 
equation, 

lmT(E,k,k') = -^ / ^T*(E,k,q E )T(E,q E ,k'), 

(21) 

where \c[ E \ = li 1 \ / mE, the integration is performed over 
the direction of this vector, and the complex conjugate 
amplitude T*(E, k,k') obeys Eq. dT2j> with -iO in the 
denominator of the integral term. This relation results 
in the unitarity condition for the scattering matrix (opti- 
cal theorem), and its validity in second order of the per- 
turbation theory is demonstrated in Appendix D. The 
analog of Eq. (|2"Tj) for partial wave scattering amplitudes 
T m (E, k, k') with azimuthal (magnetic) quantum number 
to follows from (1211) after integrating over the directions 
of k and k' with the proper angular harmonic. As an 
example, for the s-wave scattering channel with 



T s (E,k,k / ) = (T(E,kX)) 



(22) 



r a (£,k,k') « r (1) (£,k,k') g\k-k'\l (18) we obtain 



lmT s (E, k, k') = -—T* S (E, k, q E )T s (E, q E , k'). 

On the mass shell, k = k' = q E = h~ 1 \/mE, the above 
relation reads 



imr 6 .(fc) = -^|r s (fc)| 2 . 

This implies that up to the second order one has 



Imr.(fc) 

where 

r a (fc)«rW(fc) = - 

2nh 2 



Ah 2 



[Rer s (fc)] 2 , 



(23) 



(24) 



27Tfi 2 



gkl [L_i(2fc/)-Ii(2fe0] 



gkl- ( 1 - -kl 
m it V 2 



(25) 



is just angular average of Eq. (fl"4)) . In Eq. (f25|). L„(z) 
and l n (z) are the modified Struve and Bessel functions, 
respectively, and the Taylor expansion in powers of kl 
gives a good approximation for kl < 0.2. 



V. THE MANY-BODY PROBLEM 

Coming back to the many-body problem, we note that 
the amplitude of the interlayer scattering in all three 
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regimes is negative in the s-wave channel (in the regime c 
this requires E ~ Ep 3> Ej,, which is realistic in the limit 
g <C 1)- This means that at sufficiently low tempera- 
tures, the bilayer fermionic dipolar system undergoes a 
BCS pairing transition into a superfluid state with in- 
terlayer s-wave Cooper pairs, characterized by an order 
parameter A(p) ~ ^_(p)^>+(— p)^ with «/> a (p) being 
the field operator in the momentum space, which is in- 
dependent of the azimuthal angle <p, A(p) = A(p). 



A. BCS approach to pairing 

The critical temperature T c of this transition is calcu- 
lated from the linearized gap equation. In the simplest 
BCS approach, which does not take into account many- 
body effects (see below), this equation for the considered 
system is (in what follows we will use the wave vector k 
instead of the momentum p) 



A(k) 



dk' ~ tenhfc^ 
2^2d(k - k ) — A(k ), 



(27T) 



26' 



where = h 2 k 2 /2m 



(26) 

K 2 (k 2 - k'p)/2m and 
V2£>(k — k') is given explicitly by Eq. ([8j. In the regime 
^fI ^ 1, this equation can be solved directly. For a di- 
lute gas (kpl < 1) , however, the gap equation (l26l) with 
the bare interparticle interaction V2_o(k) is not conve- 
nient because it mixes many-body physics (BCS pair- 
ing) with the two-body one (scattering). In a dilute 
gas, they are well-separated in momentum space: the 
pairing originates from the momenta of the order of the 
Fermi momenta, p ~ pf, while the two-particle scat- 
tering is related to high momenta p ~ h/l 3> pf that 
correspond to short interparticle distances, at which the 
presence of other particles is irrelevant and physics is de- 
scribed by the two-particle Schrodinger equation. For 
the pairing problem, the two-body physics can be taken 
into account by expressing the bare interparticle inter- 
action V2d(p — p') in terms of the scattering amplitude 
T(E, k, k') using Eq. (TT2"1) . This results in the renormal- 
ized (linearized) gap equation 



A(k) 



dk' 

(2^)2 
1 



r(2/x,k,k') 



tanh(&-/2T c ) 



26' 



2fi- h 2 k' 2 /m + i0 



A(k'), 



(27) 



where we choose E — 2/i for convenience. The contri- 
bution to the integral in this equation comes only from 
momenta p = hk' ~ pp and, therefore, the form (|27p of 
the gap equation is more suitable to describe the BCS 
pairing in a many-body dilute system. 

In the regime of weak coupling characterized by a small 
parameter A — vpT <C 1, where vp = m/2irH 2 is the 
density of state on the Fermi surface (k = k! = kp), 
this equation can be solved by expanding in powers of 
A (see below) or numerically. However, the linearized 



BCS gap equation (|27j) can only be used for the calcula- 
tion of the leading contribution to the critical tempera- 
ture, corresponding to the terms ~ A -1 in the exponent. 
As was shown by Gor'kov and Melik-Barkhudarov [25| . 
the terms of order unity in the exponent affecting the 
preexponential factor in the expression for the critical 
temperature, is determined by the next-to-leading order 
terms, which depend on many-body effects. In the con- 
sidered fermionic dipolar system, these effects result in 
the appearance of the effective mass to* and the effec- 
tive interparticle interaction. The latter corresponds to 
the interactions between particles in a many-body system 
through the polarization of the medium - virtual creation 
of particle-hole pairs. The BCS pairing with the account 
of the many-body effects can be viewed as a pairing of 
quasiparticles of mass to* interacting with the effective 
interaction V c g. Note that, in contrast to the Fermi gas 
with a short-range interaction, in which the difference 
between the bare m and the effective to* masses (or, in 
other words, between particles and quasiparticles) are of 
the second order in A, in the dipolar system this difference 
is typically of the first order in A due to the momentum 
dependence of the dipolar interactions. 



B. The role of many-body effects 

Qualitatively the role of the many-body effects in the 
gap equation can be understood as follows. After per- 
forming the integration over momenta, the gap equation 
can be qualitatively written as 



1 = vpVea 



Inii 



C 



(28) 



where v* F — m*/2irH 2 is the density of states of quasi- 
particles with the effective mass m*, and we replace the 
scattering amplitude T with some effective interaction 
V^ff = r + SV with SV being the many-body contribution 
to the interparticle interaction. Note that the (large) log- 
arithm In fi/T c results from the integration over momenta 
near the Fermi surface, whereas the momenta far from 
the Fermi surface contribute to the constant C ~ 1. We 
can now expand v F V e n in powers of A up to the second or- 
der term, u F V e s — A + aA 2 , where the first term results 
from the direct interparticle interaction and the many- 
body effects (the difference between to* and to together 
with SV) contribute to the second term. In solving Eq. 
(|2"8f iteratively, we notice that Alnyu/T c ~ 1 and, there- 
fore, the terms aA 2 In fi/T c and AC are of the same order. 
As a result, both terms have to be taken into account for 
the calculation of the critical temperature. It is easy to 
see that they contribute to the preexponential factor in 
the expression for the critical temperature. These con- 
tributions are usually called Gorkov-Melik-Barkhudarov 
(GM) corrections [25|]. It is important to notice that the 
many-body effects appearing in Eq. (|28[) only in be com- 
bination with the logarithm In fi/T c that originates from 
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'(1) 



V- 

— V 



+ 



a 



a 



a 



i V 

a a 



+ 



CC (2 



FIG. 3: The first order diagrams for the fermionic self-energy. 
Dashed lines correspond to the dipole-dipole interactions. 
The first two diagrams contain only the intralayer interac- 
tion, while the last one describes the effects of the interlayer 
coupling. 



momenta near the Fermi surface. Therefore, it is suffi- 
cient for our purposes to consider the many-body contri- 
butions only at the Fermi surface, and the renormalized 
linearized gap equation with the account of the many- 
body effects reads 



A(k) 



TO* 

m 



dk! 



T(2 M ,k,k') 



tanh(£ fc '/2T c ) 



26' 



I 



iO 



A(k') 



(29) 



2/i — h 2 k' 2 /m 

-^^(k,kQ tanh ^ /2Tc) A(kO 



k'<Aki 



(27T) 



2^ 



where we introduce an upper cutoff Akp with A ~ 1 for 
the purpose of the convergency at large momenta. As 
discussed above, the exact value of A is not important 
because the large momenta contribution to this integral 
has to be neglected. 



It is easy to see that H^\u,p) is frequency independent, 
S« (w,p) — S a (p), and the momentum dependence re- 
sults only from the first diagram that corresponds to the 
exchange intralayer interaction (the momentum indepen- 
dent part of the self-energy leads to unessential change 
of the chemical potential) . The analytical expression for 
T,a\p) reads (p = hk) 



dq 



(2tt)2 
V 2D (0)n 

W) 



V ++ (k-q)-V ++ (0) N(q) 



V ++ (k-q)-V ++ (0) N(q) 



where N(q) = 9(kp — q) is the Fermi-Dirac distribution 
for zero temperature (the usage of n q at zero temperature 
is justified by the exponential smallness of the critical 

temperature T c ) and V-\ (k) is the Fourier transform of 

V++(p), and straightforward calculations with the usage 
of Eqs. (0 and ^ gives (see, for example, Ref. [27| ) 

to 4 4 

— = 1 - —a d k F = 1 - —gkpl. (30) 
to o7t 6ir 

It is easy to see that higher order contribution will 
introduce small parameters gkpl$ or gkpl and, therefore, 
can be neglected. 



Effective inter-particle interaction 



1. Effective mass 

The contribution to the effective mass originates from 
the momentum and frequency dependencies of the self- 
energy E Q (w,p) of fermions (see, for example, Ref. [2t|), 
to/to* = (1 + 2to<9£/gV)(1 - dT,/duj)- l \ p=PF ^ = Q. In 
the considered case, the leading contributions to the self- 
energy are shown in Fig. [3J where the fermionic Green's 
function is 

G a {uj,p) = 1 , — — -. 

I 



Let us now discuss the many-body contributions to the 
effective interparticle interaction. We consider first the 
case when the scattering of two particles with energies of 
the order of the Fermi energy corresponds to the regime 
a (g < kpl < 1), and, hence, is well-controlled by the 
Born expansion in powers of the bare interparticle inter- 
action with the small parameter A = vpT ~ gkpl. As it 
was argued above, it is sufficient to consider only the low- 
est (second order in the interparticle interactions) many- 
body contributions to the effective interaction. These 
contributions are shown in Fig. 2J and the corresponding 
analytical expressions read: 



J ^,q+k_/2 — <sq-k_/2 

. - f - m- k-/2> ^ , ++(q _ k+/2K 

J (271T £q+k_/2 - £q-k_/2 

W> - - f A * '- ± k - /2) - ~ (q - k - /2) VM±-)V' ++ (« + 

J ( 27r ) ?q+k_/2 - 4q-k_/2 

J l 27r ) ?q+k+/2 - ?q-k+/2 
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where k± = k ± k' and we keep only momentum de- 
pendent part V, i of the intralayer potential because the 
contributions of the momentum independent part of V ++ 
in SV a , SVb, and SV C cancel each other, as it should be. 
The contribution SV a (k, k') can be calculated analyti- 
cally: for k = k' = kp we have k_ <2kp and, hence, 



SV a (k,k') = -2is F V 2D (k-)Vl + (k- 
m 2nh 2 



-,g|k-k'|Z) 2 



-(glY(k-k'r, 

III 

while the other three can be computed numerically. The 
corresponding s-wave contributions are obtained by av- 
eraging over the directions of k and k' (azimuthal angles 
tp and ip' , respectively): 



SV i = (6V i (k,V)) li 



^ ^ SVi(k, k'), i — a,b,c,d 



In the regime b (exp(— 1/g 2 ) <C kpl < g < 1) the lead- 
ing contribution to the two-particle interlayer scattering 
is given by the second order Born term, and the small pa- 
rameter of the theory characterizing the interlayer scat- 
tering is A = vfT = —g 2 /4. For the intralayer scattering, 
however, the leading contribution is still given by the first 
order Born term ~ gkpl. This is because the intralayer 
scattering occurs between identical fermions and, hence, 
the dominant contribution is the p-wave one. The sec- 
ond order Born contribution in this case is proportional 
to {gkpl) 2 ln(fcF^nHsee Ref. [28|) or {gkpl) 2 \n(gkpl) for 
lo — >• (see Ref. [29|]) and can be neglected. As a result, 
the leading contributions to the effective intcrparticle in- 
teraction will be given by the same diagram from Fig. [4) 
in which all interaction V 2 d lines that connect fermionic 
lines belonging to different layer are replaced with the 
second order Born scattering amplitude T^ 2 ). It is then 
easy to see, that the contribution to <5I^(k, k') comes from 
the diagram d and equals (k = k! = kp) 



/0 (27T) 

The contribution 5Vi, can be written in the form 
2irh 2 



SV 



m 



-(gkpl) T] a , 



where r\ a = —4 and numerical calculation of the integrals 
for SVb = 5V C and 5Vd result in 

r/6 = Vc = 1-148, Vd = 0.963. 

As a result we obtain 

2irh 2 



W=—(gk F l) 2 T,m = 



Q.741 — (gkpl) 2 . (35) 
m 



5V(k,k') « v F 



' 2TTh 2 g 2 - 


2 2nh 2 


V" 


m 4 


m 


4 _ 



= v- x \ 2 



(36) 

in the regime 

r(#,k,k') 



The interlayer scattering amplitude 
c (exp(-l/.g 2 ) < kl < g < 1) is 
(4Trh 2 /m)(\n{E b /E) + in)- 1 « (4nh 2 /m) \n-\E b /E), 
similar to the scattering amplitude for a short-range po- 
tential. The corresponding small parameter is simply 
A = 2/\n(Eb/ep). (Note the conditions Re A < and 
|A| < 1 requires ep > E b .) Arguments, similar to those 
given for the regime b, lead us to the conclusion that the 
leading many-body contribution to the effective interpar- 
ticle interaction is given by the diagram d in Fig. 21 in 
which the interaction lines are replaced with the scatter- 
ing amplitude (see analogous considerations in Ref. 30] ) : 



SV 



iCC = ± 



a 

+ 



b 



+ X 
— / \ 
-9 — < — >— *- 



SV(k,k') 



2-kK 2 



m In 2 (E b /e F ) 



v- F X \ 2 . 



(37) 



Note that the leading many-body contribution to the 
effective interparticle interaction in the regimes b and c 
can be written as 



5V{k,k') = v F x \ 2 



(38) 



which is independent on the directions of k and k' and, 
hence, coincide with its s-wave component, SV — Vp 1 )? . 



FIG. 4: The second-order contributions to the effective in- 
terlayer interaction. Solid lines correspond to particles from 
different layers (labeled by + and — ) and the dashed lines cor- 
respond to dipole-dipole interactions. Note that the diagram 
d contains only the interlayer interaction, while the diagrams 
a, b, and c have both the inter- and intralayer interactions. 



VI. CRITICAL TEMPERATURE IN THE 
DILUTE LIMIT 



We now proceed with the solution of the gap equa- 
tion ((29)) in the dilute limit kl < 1. As we have already 
pointed out, the order parameter has the s-wave symme- 
try, A(k) = A(&), and, therefore, it is convenient to work 
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with the gap equation projected to the s-wave channel: 
k'dk' 



A(k) 



-o 2ir 
tanh(&//2T c ) 



T s (2fi,k,k') 



1 



AA-j 



26' 2fi-h 2 k' 2 /m 
k'dk'tanh(&/2T c ) 



iO 



2tt 



-SV- 



26' 



A(fc') 



A(fc') 



(39) 



where the first term comes from the part (a). 

It is easy to see that the first term is larger than the 
second and the third ones by a factor In [2 exp(7)w/7rT c ], 
and, therefore, the last two terms contribute only to the 
preexponential factor in the expression for the critical 
temperature. In order to solve Eq. (|42|) . we choose oj 
such that 



where r s (2/i, k, k') is the s-wave component of the vertex 
function for the interlayer scattering, see Eq. (1221) . 

Note, that following our previous discussion, the com- 
bination (m*/m)T s (2fj l) k, k') in Eq. (|39|) has to be cal- 
culated to second order in the small parameter, and the 
second order term has to be taken at the Fermi surface 
(k = k' = hp), similar to the SV contribution. All these 
second order terms can be treated perturbatively. 



A. BCS approach 

In the first order in the small parameter Eq. (|3"9")) cor- 
responds to the BCS gap equation 

l/dk '-T s (2^k,k>) (40) 



In ^R(0, -n)A(-n)+ [ d? In 

u J-u 



—[R( ,z>)A(a=o 

(43) 



and, putting £ = in (|4"2"j) . we obtain the following equa- 
tion to finding the critical temperature: 



A(0) 



, 2 expMui 7r 

In i — 

ttT c 2 



i?(0,0)A(0). (44) 



It follows from this equation that 



1 



A(fc)=- / 
Jo 



2 exjpMuj 7r 

In l— = ; r , 

ttT c 2 #(0,0)' 

and, therefore, after using Eq. 



I? cs = 2exp(7) h;exp 



J?' (0,0) 



(45) 



(46) 



2tt 

tanh(^/2T c ) 



A(k'). 



26' 2fi- h 2 k' 2 /m + i0 

In order to solve this equation, we rewrite it in the form 

"tanh(&72T c ) 1 



where B! is the real part of R, R' = KeR, such that 
R = R' + iR", and R" m -(vr/2)(i?') 2 , according to Eq. 

The value of the energy u can be obtained from Eq. 
(El: 



A(0 



26' 



26' - *o 



A(6 



(41) 

where = ^ F r s (2/x, k^, k'^,) with £ = h 2 k 2 /2m-(i, 

fcfr = Tr 1 ^/2ni{S > + /i), and &£, = h- 1 y j2m{£' + //). 

We then introduce a characteristic energy w, which is 
of the order of the Fermi energy and, on the other hand, 
is much larger than the critical temperature, u 3> T e , and 
divide the integral over £' in Eq. (j4"Tj) into three parts: 
(a) the integration of i?(£, 0)A(0) from —uj to u, (b) the 
integration of R(£, £)A{£) - R(£, 0)A(0) from -u> to u, 
and (c) the integration of i?(£, £')A(£') from —fi to — cj 
and from w to oo. In the part (a) we use the asymptotic 
formula 



'lncj = In fi 



1 



J?(Q,-M)A(-m) 

i?(0,0)A(0) i?(0,0)A(0) 



dZ']n\Z>\ — [R(0,t;')A(a 



(47) 



Substituting of Eqs. (gSJ) and O into Eq. (gSJ) results 
then in the equation for the order parameter 



/ tanh(6'/2T c 



26 



In 



2 exp(7)w 

7TT,, 



#(0,0) 



dg In 



d 



R(0,Z')-R(Z,?) 



A(0 , (48) 



while in parts (b) and (c) we replace tanh(6'/2T c ) by 
the step function (omitting the unimportant contribution 
from a narrow interval |£'| < T c -C w) and integrate by 
parts. Eq. (|4ip then takes the form 

2 exp(7)w 



in which the second term is proportional to the small pa- 
rameter and, hence, can be considered as a perturbation. 
Therefore, to leading order in the small parameter, the 
solution of Eq. (|4"5)l reads 



A(0 



In- 



i?(e,o)A(o) 



A(0 « #Sa(0). 



(49) 



In^,-/i)A(-/i) 



d£' In 



7 [i?(£,£')A(0], (42) 



i?(0,0) 

Substituting this expression into Eq. (|47p. we obtain 

(50) 
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where we can replace all functions # with their real parts 
#', see Eq. This expression, together with Eq. 

(|46|) . provide the answer for the critical temperature in 
the BCS approach. 



B. Critical temperature in the many-body system 



Following our previous discussion, one should take into 
account only those many-body contribution that appear 
in combination with the large logarithm ln(ep/T c ) orig- 
inating from the momenta close to the Fermi momenta. 
Therefore, in view of the many-body contributions, Eq. 
(|4"4")l can be written as 



A(0) 



TO* [ 2e 7 W 7T 

In i — 

to [ 7rT c 2 

ln^^A(O) 
irl r 



—R(0,0) + v F 5V 
m 



#(0,0) A(0) 



2e 1 u> 7r 

In i — 

nT r 2 



A(0). 



Therefore, 



2e 1 u> it 
In — — I— = 



ttT c 2 m*R(0,0)/m + v P SV 
1 

m*#'(0, 0)/to + iR"(0, 0) + v F W 

i (ff - 1) #'(o, o) + ^sF + i#"(o, o) 

#'(0,0) 

1 



R'(0,0) 2 



#'(0,0) #'(0,0) V to 



i /TO* \ Z/plJV .7T 

1 ) + — - I — . 



#'(0,0) 2 2' 



As a result, for the critical temperature we obtain 



T c = exp 



to*/to — 1 vp5V 



T c BCS exp 



#'(0,0) #'(0,0) 2 
to*/to — 1 i/pJV 



#'(0,0) #'(0,0) 5 



(51) 



where w and to* are given by Eqs. (ISTfl) and (|30l) . re- 
spectively. The specific expression for the many-body 
contribution to the effective interparticle interaction SV 
depends on the regime of scattering, see Eqs. (1551) . (j3"6) 
and (1371). 



We now analyze the expression (|5"Tj) for the critical 
temperature for different regimes of scattering: 



Regime a: For g < kpl < 1, we have 
a 2 

#(0,0) « !/ F r«(fe F ) - 2- {1 - 2(M) 2 [5-4 + 31n(M)]} 



15 

4 7T 

-gk F l-{\ - -M) 

7T Z 



(52) 



^- {1 - 2(M) 2 [5.4 + 31n(M)]} 4 ' 



15 ' 



TO* 



- 1 ~ —z-gkpl 

6n 



v F 8V ks -0.741(.gM) 2 , 



(53) 
(54) 



where we expand v F Ts (kp) in the expression for 
#'(0,0) up to the second order in powers of kpl, and 
keep only those terms that give contributions up to or- 
der unity in the expression for the critical temperature. 
For the calculation of u>, see Eq. (|50l) . it is sufficient to 

take #(0,0) in the form #(0,0) = v F (^(k - k')\ 

The resulting integration can be performed in the same 
way as for the integral li from Appendix D, and we ob- 
tain 



u) = fi exp 



-0.697 



(!)' 



0.651/x. 



With the help of Eqs. (j52")) -([5~4" |) we can write (within the 
accepted accuracy) 



#'(0,0) 



Up 



2 4^3 

r + is" 



-I -1 



r«(fc F 

(J) 2 [5.4 + 31n(M)] 



1 /7T\2 
2 



gkpl-(l-^kpl) + ^- 

7T I 4 

1.17- 0.925 ln(M), 



15 



(55) 



(56) 



to*/to — 1 1 
#'(0,0) ~ 3' 



[#'(0,0)] 



-0.741 



U 



-0.457. 



From Eq. (|51[) we now obtain the final expression for 
the critical temperature for the BCS pairing in the regime 
a of interparticle scattering 



2e 7 

T c a = 0.651/xexp 

7T 



"f rW(fc F 

1.17- 0.925 ln(M) 



4 15 



-1 



0.457 



0.259/z(M) 



-0.925 



(57) 



exp ■ 



gkpl 



kpl 



r + is" 
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where for kpl < 0.2 one has vpT^' {kp) « gkpl^(l — 
^kpl). The comparison of this results with the one ob- 
tained in Ref. [12| in the BCS approach shows that the 
many-body effects result in a larger (by a factor of two) 
numerical prefactor. This is the effect of the compe- 
tition of decreasing of the critical temperature because 
of the smaller effective mass and increasing T c because 
of the attractive many-body contribution to the inter- 
particle interaction. In addition, Eq. (|57|) contains an 
extra term 4g 3 /15 in the denominator in the exponent, 
originating from the third order contribution in the par- 
ticle scattering amplitude. This term is smaller than 
the other two. However, one needs a stronger condi- 
tion, namely g < (kpl) 3 / 2 , to neglect this term. This 
is because being expanded, this it results in the contri- 
bution ~ g 3 /(kpl) 3 = (g/kpl) 2 (kpl)~ 1 , which is small 
only under the stronger condition. Note that this term 
also leads to the higher critical temperature the in the 
BCS approach. 

The order parameter, Eq. (|49|) . in this regime has the 
form (k' = kp) 



A{k) ~ 5 (|k-k'|/exp(-|k-k'|0) 



4 



.9<|k 
2 



k'|z(i-|k-k'|z)>, 

Akkp 



4 



(k + k F )l E[ 



(k + k F f 



where E(z) is the complete elliptic integral and we as- 
sume kl, kpl < 0.2 to ensure the reasonable accuracy of 
the truncated expansion. 

Regime b: In the regime b, kpl < g < 1, we have 

^ 2 a ~ 



#(0,0) 



9" 4 7T . ; 

— : gkpl — (1 kpl) 

4 J F 7T 2 * ' 15 



9_ 

32 
2 



4 r 



ln(4:h 2 /mnl 2 ) 



7 



-2 7 



m 



4 7T 

i / / \ -gk F l-(l - -M), 

4 

3n 



(58) 

(59) 
(60) 



The leading order contribution in #(0,0) is momentum 
and energy independent and, therefore, we have io = [±. 
From Eqs. (|58" 1) -(|60" 1) we obtain 



1 



#'(0,0) 



— + gkpl-(l - -kpl) + — 



4 7T 

4ti 2 



15 



hi 



m*/m — 1 14 A^Z 



#(0,0) 
[#(0,0)] 5 



3tt ,g 
1. 



< 1, 



^-2 7 
2 ' 



(61) 



Note that the many-body contribution to the effective 
mass is negligible because the leading term in the scat- 
tering amplitude is momentum and energy independent. 
From Eq. (|5ip we then obtain 



rp 2e^ 

T cb = ^iexp ■ 

1 

+ 2 



In 



9 2 

T + 
\to/zZ 5 



gk F l-(l 

7T 



-2 7 



-k F l) 



- 1 



15 



(62) 



4 e 3/4 



/ifi 2 

to/ 2 



gkpl ( 2fci?/ 



4 5 3 
15" 



2.7 



M 2 

^2 ex P 

Till 



g 2 /A + Agkpl/ir - 2gk 2 F l 2 + 4 ff 3 /15 



Furthermore, we note that the exponent in this expres- 
sion coincides with that in Eq. (j57| and, following the 
same arguments as in Eq. (|57[) , we keep some higher 
terms in the denominator in the exponent. 

Regime c: Finally, for exp(— 1/g 2 ) < kpl <C g < 1, 



#(0,0) 
1 



rn 



HEb/lA ' 
4 

gkpl <C 1, 



3tt 



where is given by Eq. (jllj) . Therefore, taking into 
account that, similar to the regime b, lj = fx, we obtain 



^exp^ -]n(E b /fj,) 



2exp(7 — 1) 



(63) 



This expression is completely analogous to the critical 
temperature in a two-component 2D Fermi gas with a 
short-range interparticle interaction (see Ref. [30]), as it 
should be in this regime. 

Note that within the accepted accuracy, both expres- 
sions (1521 and for the critical temperature in the 
regimes b and c can be written in the form 



rp 2e^ 

J c ,bc = yuexp 

7re 



0.42^ exp ■ 



2 gkpl 



kpl 
2/V 



(64) 



HEb/fi) 



4 7T 

gk F l-(l - -kpl) 

7T Z 



In both these regimes, the order parameter is to the lead- 
ing order momentum independent, A(fc) ~ const. 

Eqs. (|57|) and (|64|) provides the answer for the criti- 
cal temperature of the BCS transition in a dilute bilayer 
dipolar gas. The corresponding values of T c calculated 
according to these formulae are small, T c < 10 _3 -^10~ 2 /i, 
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because the exponent in Eqs. (|5T[) and (1641) contains the 
inverse of the product (or square) of the small param- 
eters of the problem. For example, even for kpl = 0.5 
and g — 0.45 < kpl one has T c ss 10 _2 /u. This makes 
an experimental realization of the superfluid state very 
challenging. The situation is more promising in the dense 
case. 



VII. CRITICAL TEMPERATURE IN THE 
DENSE LIMIT 

We assume now that kpl > 1 and j « 1 such that 
gkpl — kpdd < l(and kplo <C 1). The condition gkpl < 
1 ensures the validity of the perturbative expansion in 
powers of the interlayer interaction, although the mean 
interparticle interaction is comparable or larger that the 
range of the potential I. For this reason there is no need 
to renormalize the gap equation: two colliding particles 
are no more well-separated from the rest of the system, 
but many-particles collisions are still well-controlled by 
the small parameter ddkp = gkpl < 1. With the account 
of the many-body effects, the gap equation for the pairing 
in the s-wave channel reads 



m Jo 



2tt 



26' 



(65) 

where V e e^ s (k, k') is the effective interparticle interaction 
in the s-wave channel, 



V eS , s (k,k') = (V 2 D(k-k') + 5V(k,k' 

= Ti 1 Hk,k') + (SV(k,k')), 



rM(k,k') = (v 2D (k-k') 



Here 



and, as before, all many-body corrections have to be 
taken at the Fermi surface (assuming T c <C (J-)- Using 
the same arguments as in the previous Section, we can 
argue that the critical temperature T c is related to the 
critical temperature in the BCS-approach T® cs (with no 
many-body contributions) as 



T c = T C BCS exp 



m*/m — 1 vpSV 



i?'(0,0) [R'(0,Q)] a 



(66) 



where i?'(0, 0) = vpT^ {k F ,k F ) and T BCS is determined 
by the BCS gap equation 

r«M/v 2D (k _ k0 \ 

J Ztt \ I <p, v > 

tanh(6/2T BCS ) 



00 k'dk' 



2n 



r«(fc,fc 



A ( fc ') 

2<U" 

, tanh(^/2T c BCS ) 



26 



A(fc'). 
(67) 



The solution of this equation is actually given by Eqs. 
(|46l) and ([50j) : although we are dealing with the non- 
renormalized gap equation, the renormalization can still 
be performed as a formal trick. [Note, that within the ac- 
cepted accuracy, R'{0, 0) in Eq. (|46|) has to be calculated 
up to the second Born approximation, while only up to 
the first order in Eq. (|50]).] Nevertheless, we give here 
an alternative solution of the gap equation that follows 
the lines of Ref. 13 ill and avoid the renormalization. 



A. BCS approach 

We rewrite Eq. (f6"T| in the form 
f ,tanh(e72T BCS ^ 



/OO 
-u 



2 £' 

where £ = h 2 (k 2 — k F )/2m, 



R(t,?)W), (68) 



R(£, £') = u F r™ (k, k') = gl f ^ e -Vfe 2 +fc' 2 -2fcfc'cos V 

JO 7T 

x ^/k 2 + k' 2 -2kk'cosip, (69) 

and, following the method of Ref. [3l|, decompose the 
interaction function £') into a separable part and a 
remainder ?"(£,£') that vanishes when either argument is 
on the Fermi surface: 

i?(£,O = ^(0,0M0«(O + ^O (70) 

with = fl(£,0)/iZ(0,0) = R(0,£)/R(0,0) and 

r(f,Q) = r(0,f) = 0. Note that v(0) = 1 and v(£) 
decays exponentially at large momenta kl ^> 1, i.e. 
ml 2 i/h 2 > 1. Eq. (JSHJ) then takes the form: 

a(o = - r(o, oho r de tanHC ^ CS) v(OA(o 



2f 

,tanh(£'/2T BCS ) 



2? 

For £ = this equation reduces to 



r^OHO- (71) 



/OO 
-« 



,ta a h(^/27?CB) |>( 



2£' 



and, therefore, we can rewrite Eq. (|7Tj) as follows 
A(0 = ^)A(0)-/;^^^ C " 



«(OA(0) 



fi 2ie 



-r(£,£')A(£'), 



(72) 

r(£,OA(0 
(73) 



where we replace tanh(£'/2T B ) with sign(£') assuming 
that T c fj, and neglecting exponentially small contri- 
butions [the integral is converging because r(£,0) = 0]. 
In Eq. (|72[) we can now single out the large logarithmic 
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contribution that comes from momenta near the Fermi 
surface. This can be achieved by writing d£' '/£' = <f(ln£') 
and integrating by part with the following result 

A(0) = -R(0,0) |^(-/x)A(-/x)ln/i 

2 J , deinWl Nnh(e72r c BCS ) V (e')A(e') 
After performing the derivative, 



_d_ 



tanh( 



2TBCS 



)»(OA(0 



2T C BCS cosh 2 (f/2r c BCS ) 



and using the fact that 1/2T C BCS cosh 2 (£'/2T BCS ) is 
sharply peaked at £' = 0, we obtain 

f 1 2e 7 
A(0) = -i?(0, 0) |- V (- M )A(-M) In /x + In A(0) 

-5 |^e'ln|e'|sign(^)A[ w (e')A(^)] 
= - J R(0,0){ln-^A(0) 



where we replace again tanh(^'/2T ( BCS ) with sign(£'). 

The pair of equations ([74")) and ([75)) can now be solved 
iteratively because the integrals in both equations pro- 
vides small corrections when kpl > 1 (see below). In the 
leading order, we have for the order parameter<£' 

A(0««(OA(0). 

[Note that the second iteration of Eq. ([73)) 

with the explicit expression r(£, £') = £') — 

i?(C,0) J R(0,^)/ i ?(0,0) can be rewritten in the form of 
Eq. (|I8|).] From Eq. we then obtain 



^BCS 



2e^ 



-/x exp 



exp 



R(0,0) 



(75) 



for the critical temperature in the BCS approach. 

To establish the connection with the approach from 
the previous Section, we notice that 



'\21 



1 PVi l^'l d r r^i x T ^'1 l^'l d r 



The last term in the right-hand side can now be identified 
with the contribution to lnw, see Eq. (1501) . while the first 
one with the second order Born contribution to i?'(0, 0). 



Before going further, let us discuss the conditions for 
the iterative approach to the system of Eqs. ([74"]) and 
([75)1 to be legitimate. In Eq. (1751) . this requires that 
the second term is small and, hence, A(£)/A(0) « 
For kp I ~ 1, one can see that the relative contribution 
of the second term is of the order of gkpl. Therefore, 
the iterative scheme with the result A(£) ~ v(£)A(0) 
is legitimate for gkpl — a^l < 1. For a dilute system 
with kpl <C 1, the situation is more subtle. In this case, 
the contribution from £' < /x in the second term gives 
the relative contribution of the order of gkpl, while the 
regime /x < £' < h 2 /ml 2 results in the relative contribu- 
tion ~ g/kpl. Therefore, for g < kpl both contributions 
are small and the iterative procedure is legitimate. How- 
ever, for the opposite case g > kpl, the contribution from 
the second region is large and the iterative scheme breaks 
down. The reason for this is the large value of r(£, £') for 

3> /i. The proper iterative procedure in this case can 
be developed on the basis of the renormalized gap equa- 
tion. 



B. Critical temperature in the many-body system 

The calculation of the many-body contributions to the 
critical temperature can be performed in the same way 
as in the dilute regime because, although kpl > 1, we 
assume that gkpl — a^kp < 1, and, hence, the expansion 
in powers of interparticle interaction is still valid. The 
corresponding contributions to the self-energy (effective 
mass) and interparticle interaction are shown in Figs. [3] 
and0J and the analytic expressions are given by Eqs. ([50)) 
and ([3i"|) - ([34]) . respectively. Note that 514 (k, k') can be 
calculated analytically, 

514 (k, k') - -2v F V 2D (k - f)V| + (k - k') 
Airfr 2 

« ( ff |k-k'|0 2 exp(-|k-k'|0, 

m 

[as before, we keep only the momentum dependent part 
of V ++ (k— k')] together with its angular average for k — 
k! = kp, 



u F SV a = - 8(gk F l) 2 |l (2M) - g 0*1 (2; kpl 2 ) 
-^(M)iF 2 (2;||^ 2 ) 



where 1q{z) is the modified Bessel function and 
n F m [a\, . . . , a n ; b\, . ■ ■ , b m ; z) is the hypergeometric func- 
tion, while the other three contributions require numeri- 
cal integrations. We write these contributions as 

upSVi = (gk F l) 2 fi(kpl), i = a,b,c,d, 

where the functions fi(x) are shown in Fig. [5] The over- 
all angular averaged contribution to the effective interac- 
tion then reads 

vpSV = (gkpl) 2 f (kpl), 
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fi(x) 



\ \ 
\ \ 

\ \ 



•"'2 



4 X 5 



FIG. 5: The functions f a (x)/4, f b (x) = f c (x), fd(x), and f(x) 
(dotted, dash-dotted, dashed, and solid lines, respectively). 




FIG. 6: The function -y{x). 



where the function f(x) = f a {x) + 2ft,(x) + f c (x) is also 
shown in Fig. [5] (solid line). 

After writing R(0, 0) in the form 



#(0,0) = — gk F lMk F l), 

7T 



where 



7 {x) = \ £ dtp Bm{<p)e-' sin ^ = \ [L_ a (2x) - h (2a:)] , 
see Eq. (|25p . is shown in Fig. [5J we obtain 



1 (m* \ vpSV 



R(0,0) V m V [i?(0,0)] 2 

1 /7TX2 /(fc F 



(76) 



3 7 (M) V4/ 7 (M) 2 ' 
Finally, after combining together Eqs. (|75|) . (|6"6")l and 
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/ / r(0.7,M) 
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1.5 



FIG. 7: The function r(g,x) for g = 0.7 (solid line), g = 0.8 
(short-dashed line), and g = 0.9 (long-dashed line). 



([TBI, we find 



T c = ^iexp 

ir 



exp 
where 



1 n(fc F Z) 1 
167(M) 2 ~ 3t(M) 
1 



V4 



4gk F l j(k F l ) 



7(M) 2 . 
(77) 



1 f°° <f 

= - / rfsln|l-s 2 |sign(s-l)— [F(s,a;) 2 ] 

2 Jo 



and 

v( s ,z) 



2s cos (pe 



-X-y/l + S 2 — 2s COS 7p 



The expression ((77)) is appropriate for k F l > 1. Based on 
the discussion after Eq. (|75)l . we can write the expres- 
sion for the critical temperature that will interpolate the 
behavior for kpl > 1 and k F l < 1: 



T c = ^exp 

7T 



exp 
2eT/x 



1 



3 7 (M) 

7T 



7(M) 2 
l 



4 5 M 7 (M) 1 - (4/tt) 5 M7(MMM) 
r(ff,M). (78) 



The dependence of the function r{g,k F l) on k F l for 
several values of g is shown in Fig. [7J We see that the 
critical temperature decreases very rapidly for k F l > 1 
due to the fast decay of the scattering amplitude. The 
optimal value of k F l is around 0.5 with the critical tem- 
perature reaching values of the order of O.l/i for g ss 0.9 
that corresponds to gk F l « 0.45 < 1. 
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VIII. CONCLUDING REMARKS 
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We obtain our results for the superfluid critical tem- 
perature using the mean-field approach. However, as it 
is well-known, this approach in two dimensions is only 
applicable at zero temperature, while at finite temper- 
ature the long-range order is destroyed by phase fluc- 
tuations and, therefore, the mean-field order parameter 
is zero. In this case, the transition into the superfluid 
phase follows the Berezinskii-Kosterlitz-Thouless (BKT) 
scenario [32l l33l|. In the weak coupling limit, however, as 
it was pointed out by Miyake Q , the difference between 
the critical temperature calculated within the mean-field 
approach T c and the critical temperature of the BKT 
transition Tbkt can be estimated as T c — Tbkt ~ T 2 /pi 
and, therefore, small as compared to T c . As a result, 
our mean-field calculations provide a reliable answer for 
the critical temperature in the considered weak coupling 
regime adkp < 1. 

Let us now discuss possible physical realizations of 
the interlayer pairing. In the experiments with polar 
molecules, the values of the dipolar length ad are of 
the order of 10 2 10 4 nm: for a 40 K 87 Rb with cur- 
rently available d ss 0.3 D one has ad ~ 170 nm (with 
ad ~ 600 nm for the maximum value d ~ 0.566 D), 
and for 6 Li 133 Cs with a tunable dipole moment from 
d = 0.35 D to d — 1.3 D (in an external electric field 
~ lkV/cm) the value of ad varies from ~ 260 nm to 
ad ~ 3500 nm. For the interlayer separation I of the order 
of few hundreds nanometers, the corresponding values of 
the parameter g can be both smaller and larger than 
unity (g < 10). 

The values of the parameter kpl are also within this 
range for densities n = 10 6 -j- 10 9 cm -2 (for example, one 
has UfI = 1 for I — 500 nm and n 3 • 10 7 cm~ 2 ). Note, 
however, that the optimal values of this parameter are 
around kpl ~ 0.5 (see Fig. [7]) and, hence, the optimum 
value of the interlayer separtion is related to the den- 
sity, which, in turn, should be large enough to provide 
a substantial value for the Fermi energy. For 40 K 87 Rb 
molecules at the density n w 4 • 10 s cm~ 2 in each layer 
one has Ef ~ 100 nK and kp =. Therefore, the inter- 
layer separation I should be relatively small, I < 150 nm, 
to meet the optimual conditions. For I = 150 nm one 
then has g rs 1.1 (with current d rs 0.3 D), kpl ri 1, and 
T c O.lep R! 10 nK. Note that stricktly speaking these 
values of parameters g and kpl do not correspond to the 
weak coupling regime considered in this paper, rather 
to the intermediate regime of the BCS-BEC crossover. 
However, based on the experience with the BEC-BCS 
crossover in two-component atomic fermionic mixtures, 
in which the critical temrature confines to grow when 
approaching the crossover region from the BCS side, we 
could expect that the above value of the critical tem- 
perature provides a good estimate for the onset of the 
superfluidity in the intermediate coupling regime. 
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Appendix A: Interlayer bound state 

We present in this appendix some details of the 
calculation of the (intra-layer) bound state properties 
for small couplings jCl Our starting point is Eq. (8), 
i.e. the equation for the radial wave-function Xm z (p) 
of the bound state with binding energy Eb- Since for 
g -C 1 one has merely one (shallow) bound state of axial 
symmetry, in the following we focus on the the axial 
symmetric case, m z = 0, and derive its binding energy 
and wave-function within a series in 1/g. 

On one hand, we see from Eq. (TiTJ|) that at sufficiently 
large distances, cf. p 3> p*, we can neglect the interaction 
potential V2d(p) aud the wave- function takes the form 



Xo(p) ~ CK (y/mEbp/h), 



(Al) 



with C a constant and Kq(z) the modified Bessel function 
of the second kind and the distance p* ~ (d 2 /Eb) 1 ^ 3 — 
(glh 2 /mE b ) 1/3 » I for j » ml 2 E b /h 2 . Since > 
%/ 1 \JmEb = Pk, for Eb <C h 2 g/ml 2 , we can expand 
Eq. (jAll) for distances p <C Pk as 



Xo(p) 



\yJmEbpJ 



with 7 Ri 0.5772 the Euler constant and in particular 



p-rHxo(p)} 

dp 



In 



/ 2h \ 



Xy/mEbp, 



7 



(A2) 



On the other hand, for sufficiently small distances, cf. 
p <C p*, and weak coupling g <C 1, we can neglect the 
bound state energy, i.e. we assume Ej, <C h 2 g/ml 2 , and 
expand the wavefunction in powers of g as 



Xo(p) ~ A/" 



n=l 



with Af an overall normalization constant. Then Eq. (fTTJ|) 
gives a set of differential equations for the various terms, 



pdp P dp Xo {P) ~ h 2 g Xo {P) 



(A3) 



with Xo (fi) = an d the boundary condition Xq 1 '^) = 
5 n o. From Eq. (|A3[) we obtain the zeroth order term is 
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then a constant, i.e. X^ip) — lj while the higher order 
terms are obtain by iterated integration of Eq. (|A3I) as 



[* P2 l{pl - 2l 2 ) („_!) 

o Jo 



l pl {p 2 + l 2 f /2 ' 



The terms up to forth order are (with z = \/ p 2 + I 2 /I) 
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(>4 



with Li„(z) = Y^kLi z k /k n the polylogarithm function. 
We truncate the latter expansion at order g n and have 



k=0 



A^(P), 



which for p » I (cf. z ^> 1) gives 
A a Q) (p)=0, 
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FIG. 8: Binding energy (scaled), g 2 ln(ml 2 Eb/Ah 2 ), as a func- 
tion of the coupling g. Show are the numerical result E£ (solid 

line) and the analytical results E^ in n = 2, 3, 4-th order per- 
turbation theory (dotted, dot-dashed, dashed lines), as well 
as the expansion up to O(g) for n = 3, 4 (thin dot-dashed, 
thin dashed lines) from Eq. (|A6[) . Note that the numerical 
results for n — 3 and n = 4 with this scale are practically 
indistinguishable. See, however, Fig. [9] 



Comparing the latter with the asymptotic behavior 
from Eq. (|A2[) . we see that for n > 2 we can match their 
leading ~ l/log(p) behavior via the binding energy and 
thus obtain (up to order n), respectively, 



In 



In 



ml 2 E [2) 
4h 2 e- 2 ~< 



(3) 



In 



ml 2 E, 



Ah 2 e- 2 ~t 



ml 2 Ef ] 
4h 2 e- 2 ^ 



9 9 
8/ff 2 - 8/g 



5, (A6a) 
5 321n(2)-ll 



1+9/15 

128 1253 
+ 15g ~ 225 



15/g- 
0(g), 



1 



(A6b) 



8/g 2 - 8/g + 5 - ll.g/15 - 311g 2 /240 
l + <?/15-.g 2 /240 

32g(ff+l)ln(2)/15-7r 2 g 2 /48 

l + g/15-# 2 /240 
8 128 2521 ^. . , 4 „ . 

? + T^-l5b- + 0(5) - (A6c) 



We see that the truncation at order n yields the correct 
expansion for the logarithm of the energy up to 0(g n ^ 1 ), 
and thereby obtain from the latter obtain 



Eh 



Ah 2 

ml 



2 ex P 



8 128 2521 „ 
7 2+ Wg-^-^ + ° { *\ 



Concluding, we compare the analytic results for the 
binding energy of Eq. (|A6|) , , with the result obtain 
by numerical integration of the Schrodinger Eq. (|10[) . E£, 
which are shown in Fig. [5] and Fig. [HI 
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FIG. 9: Deviation of the analytical results from the numerical 
result for the binding energy, _B^ n '/_B°, as a function of the 
coupling g for n = 2,3,4 (dotted, dot-dashed, dashed lines). 
Also show are the expansion up to O(g) for n — 3,4 (thin 
dot-dashed, thin dashed lines) from Eq. (|A6[l , 



Appendix B: Low energy scattering 

In the following we discuss the connection of the 
bound-state and in particular its binding energy E b with 
the scattering amplitude at (very) low energy for weak 
coupling g <C 1. Our starting is the scattering wave func- 
tion tp^ (p) a t energy E = Hk 2 /m of Section 3. 

On one hand, for p 3> it takes the asymptotic form 



^i +) (p) ~ exp(ikp) 



ifk 



with fk the scattering amplitude and Hq 1 \z) is the Han- 
kel function. In the regime kp+ <C 1 we can further ex- 
pand ip^' (p) for p <C l/k as 



-7-y 



(Bl) 



On the other hand, for p -C p* and weak coupling 
g -C 1 we can neglect the energy, i.e. take E — > 0, and 
expand the wave-function as a power series in g as 



E 

.71=0 



9 Xo (p) 



(B2) 



where M is a normalization constant. The individual 
terms, Xo^[p)i are normalized as Xo (0) = ^n,o an d are 
derived in Appendix lAl up to order n = 4, see Eq. (|A6|) . 
For p 3> Z (but still p <C p+) the wave- function approaches 
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log(p/0 



960 



which has the form of Eq. (|B1I) and thus we can match 
the two asymptotic forms. This is conveniently done in 
terms of the log-derivative of the wave-function as 



p§- p H4 +) (p)} 
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+ 7 + 
A( 5 ) 



2vr 
Ik 



in 

y 
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where A(<?) is obtained as in Eq. (|A5 
-8/.9 2 + 8/.g - 5 + O(g) 



(B4) 
(B5) 



as 



A( 5 ) 



1 - .g/15 + .g 2 /240 + O(g 3 ) 
128 2521 
g" \hg 450 



ml 2 E b 
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which we recognize as the perturbative expansion for the 
binding energy. Matching Eq. (|B4|) and Eq. (IB5|) , we get 
the scattering amplitude explicitly as 



fk 



2n 



-in 



log(2/W) - 7 + A(.g)/2 + 2 tt/2 ]og(E b /E) + wr ' 



which recovers the universal low-energy behavior of two- 
dimensional scattering. Finally, we remark that expand- 
ing the latter expression for the scattering amplitude as 
a power series up to forth order in g we obtain 
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Appendix C: Born series for the s-wave scattering. 

In the following we derive the (s-wave) scattering am- 
plitude within a Born-expansion. We recall the relation 
of the s-wave scattering amplitude fk in terms of the ver- 
tex function r(£,k,k') from Eq. CO 



/ fc -y 2^/*(v)- — wj^ — -^E r i J ( fc )' 



with k, k' on the mass shell, i.e. k = k' = y/ mE/H 2 , the 
averaging is performed over their azimuthal angles <p and 
ip', and the contributions r^(fe) follow from the Born 
expansion of the vertex-function T(E, k, k'), cf. Eq. (|13l) . 
For convenience we introduce the s-wave potential 

V s {qi,q 2 ) = (VzdIhi -q 2 )) ¥ > 1 , ¥ > 2 = (VaD(qi - q2)} ¥ , 1 

= y ^^2D (^J q\ + ql - 2 qi q 2 cos ip 
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^h 2 l d I dtp e _ ly f 



I q\+q\—1q\qi cos if 



dl J 2vr 

which in particular for q% = (or g2 = 0) gives 



. (ci) 



V s (q,0) = V s (0,q) = V 2 u(q) = -9 
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~ ql ql, (C2) 
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FIG. 10: Born-series for the s-wave scattering amplitude. 
Shown are (a) the real and (b) imaginary part of the con- 
tributions Ts (k) of order n = 1,2,3,4 (solid, dashed, dash- 
dotted, dotted lines) as a function of momentum k. 



while for equal momenta, q\ = q 2 = q, gives 
2tt H 2 g t 



V s (q,q) = - 



-2ql[^ 1 (2ql)-I 1 (2ql)], (C3) 



with L„(z) the modified Struve function and I n (z) mod- 
ified Bessel function of the first kind. 

The first order contribution then explicitly gives 



r«(fe)Hv^(k-k')>^ 2 

2irh 2 g 



V s {k,k) 

2kl[L^(2ql)- h(2ql)], (C4) 



which vanishes for k — > 0, is negative for k > with its 
minimum of pa —0.3136 x 2Trh 2 g/m at k « 0.7131 /Z, see 
Fig. [10] (solid line), and for low- momenta, k -C 1/2, gives 



rW(*) 



2irh 2 g 



Akl 
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2(fc2) 2 + 0(fc 3 ) 
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The second order contribution to s-wave scattering is 
dq (F 2fl (k-q)V 2fl (q-k')W 



ri 2) (fc) = 



(27T) 



E -H 2 q 2 /m + i0+ 



The real part, cf. the principal value integral, in general 
has to be evaluated numerically and is shown in Fig. llOf a) 



(dashed line). We remark that it has two extrema, cf. 
pa -0.2603 x 2irh 2 g 2 /m at k pa 0.1364/2 and pa 0.0103 x 
2nh 2 g 2 /to at fc pa 2.1137/2. Moreover for small momenta, 
fc -C 1/2, the leading contribution is 
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which in particular is finite and negative for fc = 0. For 
the imaginary part we have explicitly 
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g 2 (kl) 2 2n [L_ 1 (2g2)-/ 1 (2g2)]^ 
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with a minimum of ~ —0.1544 x 2irh 2 g 2 /m at fc pa 
0.7131/2, see Fig. HUTb) (dashed line), and for small mo- 
menta, fc <C 1/2, vanishes as 



Im 



ri 2 )(fc) 



16ft 2 



3 2 (fcZ) 2 + 0(fc 3 ). 
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The third order contribution to s-wave scattering is 
convenitently splits into its real an imaginary part as 



ri 3 )(fc) = 
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(V2p(k- qi)V2p(qi - q2)V2£>(q2 - k ')) y , y 
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which are shown in Fig. QJJ] (dash-dotted lines). We notice 
that the real part has its minimum at fc = 0, where 
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where for the s-wave potential we used the representation 

/>oo 

V a {q\,q2) = 2-7T / pdpV 2D (p)J (qip)Jo(q 2 p), 
Jo 

and for the convolution of the Bessel function the relation 

r°° i 

/ dqe-* i J (q P )= === , 
Jo yV 2 + 2 2 
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and its maximum sa 0.0636 x 2nft 2 g 3 jra at k ps 0.9205//. 
While the real part of T^\k) is finite and negative for 
k — 0, we notice that its imaginary part vanishes as 



and two extrema, sa —0.1857 x 2irh 2 g 3 /m at k 
and fa 0.0059 x 2nft 2 g 3 /m at k « 1.9642//. 



0.3623/; 
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The forth order contribution to s-wave scattering is 
conveniently split into its real and imaginary part as 



r«(*) = / 



dq 1 dq 2 dq 3 (V 2 n(k - qi)F 2D (qi - qa^nfqa - qs)V 2 D(qa - k')) 
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which arc shown in Fig. [TU] (dotted lines) . We notice that 
the imaginary part of T^\k) in the limit k — > is finite, 
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f qdq V s (0,q)V s (g,0) 
J 2tt E-ft 2 q 2 /m 



2nft 2 g A TT 
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with two extrema, pa -0.2184 x 2irft 2 /m at k sa 1.0532/Z 
and sa 0.0245 x 2nft 2 /m at k sa 1.0532//. The real part 

of rf \k) diverges for k — > as (for k <C l/l) 
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O(k) 
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an has two extrema, i.e. sa 0.0937 x 2-Kft 2 g A /m at k sa 
0.5574// and sa -0.0027 x 2tt/i 2 5 4 /to at k sa 1.8438//. 
Concluding we remark that summing up the contribu- 



tions up to forth order in g, we obtain for k <C l/l, 
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which coincides with the expansion of the scattering am- 
plitude in terms of the binding energy up to forth order 
in g from Eq. (|B6[) . 



Appendix D: The s-wave on-shell scattering 
amplitude in the second Born approximation 

We present in this appendix some details of the calcu- 
lations of the s-wave on-shell scattering amplitude in the 
second Born approximation. Our starting expression is 



dq V 2Z5 (k-q)V r 2I? (q-k') 



(2tt) 2 E - q 2 ft 2 /m + i0 



where V 2 d(^ — q) is given by Eq. flSJ), k = k' = \pmEjft, 
and we have to perform averaging over the direction of k 
and k' (azimuthal angles ip and ip', respectively) 

in order to obtain the s-wave contribution. 

The calculation of the imaginary part is simple and can 
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be performed without the on-shell condition k = k' = qE- 
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= -4^ y ^^2D(k-q £ )F 2D (q B -k'), 
where qE — %/ mE/h. The s-wave contribution then is 
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On the mass shell k = k' = qE = V mE/h and under the 
condition M< 1, we have 
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and, therefore, 
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The calculation of the real part 

dq (V^(k-q B ))^(T/ 2D (q-k')V 



h 2 (k 2 — q 2 )/m 



(D2) 

where f denotes the principal value of the integral, is 
technically more involved. After introducing the new di- 
mensionless integration variable y = q/k, Eq. (|D2j) reads 
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2irh 2 



00 ji 



£ ydy 
1-2/ 2 



R 1 R 2 e-< R i +R ^ 



where e = kl <C 1, R4 = \J 1 + y 2 — 2y cos <pi , </?i = y, 
and ip 2 = v 3 '- After integrating by part we obtain 
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^(i? 1 i? 2 exp[- £ (i? 1 +i? 2 )]) vu 



(2) 

and the calculation of RefTi (k)} reduces to the calcula- 
tion of the integral 



r°° d 1 

1(e) = e 2 / dj/]n( 1 - y 2 |) — (i?ii? 2 e- 
Jo 1 dy \ 



up to the terms ~ e 2 lne. 



e(i?i+fl 2 ) 



It is convenient to split 1(e) into two parts, i.e 1(e) — 
h(e) + /2(e), where 

h(s)=e 2 r 'dyHy 2 -1)± (R 1 R 2 e-^+ R ^) , 

/ 2 (e) =£ 2 t dyHl-y 2 )^(R l R 2 e-^ R A . 

Ji dy \ 1 V i V 2 

Within the chosen accuracy, the second integral I 2 (e) we 
can be simplified 

He) « e 2 ^ dyln(l - y 2 )^ (M 2 >^ 2 

because it converges when e — > in the integrand. The 
result of averaging over angles (p\ and (p 2 can now be 
expressed in terms of complete elliptic integrals E(fc) and 
K(fc): 

h(e)=e 2 \ [ dyln(l-y 2 )^-E(k) 

1" Jo y 

[(y + l)E(fc) + (y-l)K(fc)], (D3) 

where fc = 4y/(l + y) 2 . The numerical calculation of the 
above integral gives 



h(e) ~ -0.697e 2 



(D4) 



In order to calculate the integral I 2 (e) we differential 
first with respect to y: 



h(e) = e 2 dyHy 2 -i)( \y- c ^ + y-™^ 2 



-e [(y - cos(pi)R 2 + (y - cosy 2 )i?i]} 



= ha(e) + I 2b (e), (D5) 



where we split the integral into the two contributions 
ha(s) = e 2 J™dy\n(y 2 1) (e-«*+«*) 



y — cos ip 1 y — cos ip 2 
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R\jR 2 R 2 /R\ 
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[(y - cos tp 1 )R 2 + Rx(y - cos ip 2 )}) vilfi2 

For y > 1 we can write 
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+ 0(y- 2 ) 1 (D6) 



(y - cosip 1 )R 2 + Ri(y - cosy 2 ) « 2y 2 - 2y cos 



sin ipi + sin ip 2 + 4 cos ipi cos y 2 



+ 0(y' 1 ). (D7) 
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Note that the integral /26(e) is already proportional to e 3 
and, hence, the contribution of finite y (1 ~ y -C e _1 ) will 
be beyond the necessary accuracy. Therefore, we should 
consider only the contribution of large y (y ~ e -1 ). In 
this case we can replace ln(y 2 — 1) with 21ny and use 
Eqs. flD6j) and (fD7j) . In the exponent exp[-e(i?i + i? 2 )] 
we keep only 2ey from the expansion for e(R\ + R2) to 
ensure convergency, while expand in e(cos<^i + cos (^2) 
and (£/2y)(sin 2 </?i + sin 2 1^2) to the second and the first 
order, respectively, 

e -e(R 1+ R 2 ) _ e -2ey ^ + £ ^ os ^ + CQS ^ 

e 2 e 
+ — -(cosyji + cos<p 2 ) 2 - 7— (sin 2 + sin 2 ip 2 ) + ■■•]■ 
2 2y 

It is easy to see that higher order terms in the above 
expansion, as well as the omitted terms in Eq. (|D7|) 
result in terms ~ e 3 lne or smaller. The integrations 
over y and the angles tp± and (p 2 are then straightforward 
and give 



/26(e) «--+7+ln(2e)+£ 2 
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where 7 ~ 0.5772 is the Euler constant. 

The integral /20(e) can be rewritten in the form 
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In the second integral, /202(e), we can write [cf. Eq. (|D6I) ] 
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The calculations are now similar to those that lead us to 
Eq. (|D8I) . and we obtain 



:~ J dyln(y z -l) 
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R1/R2 R2/R1 



-[l-37 + 31n(2e) 



2 

In order to calculate the first integral, /201(e), in 



Eq. (|D9|) . we notice that for y 3> 1 
y — cos ipi y — cos ^2 
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and rewrite /201(e) as follows 
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The second and the third integrals in the right-hand-site 
of this equation converge at y ~ 1 and, therefore, we can 
replace e~ 2£y with unity. As a result we obtain 
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The integrations over y and angles ipi and </?2 are then 
straightforward and we obtain 



y - cos tpi y - cos <p 2 



dy\n(y 2 - l)e~ 2ey 
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where again fe = 4y/(l + y) 2 . Numerical evaluation of 
the remaining integral gives 
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and, hence, 

We) ~ 1 - 7 - ln(2e) + 2e 2 [ 7 + ln(2e) + 0.0192] 

+ ^£ 2 [1 - 57- 51n(2e)] (D13) 
= 1 - 7 - ln(2e) + e 2 [2 7 + 21n(2e) + 0.03834]. 
Combining together Eqs. (|D4j) . (|D13|) , and (|D8]l . we get 
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